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We report that the triploid apogamous species Dryopteris hondoensis Koidz. shares its sequences with 
the D. erythrosora complex (D. erythrosora, D. caudipinna and D. koidztimiana) and D. chinensis at two 
nuclear markers (PgiC and GapCp). This is the first report of reticulate evolution between subg. Erythr- 
ovariae sect. Erythrovariae, to which D. erythrosora, D. caudipinna and D. koidzumiana and D. hon¬ 
doensis belong, and sect. Aemulae of subg. Dryopteris, to which D. chinensis belongs. Apogamous spe¬ 
cies of Dryopteris may hybridize with distantly related species congeners. 
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Apogamy, or agamospory in ferns, is a type of 
asexual reproduction in which unreduced spores 
are formed. The resultant gametophytes produce 
sporophytes of the next generation without fertil¬ 
ization (Manton 1950). Apogamous reproduction 
is common in ferns. Approximately 10% of all 
species of ferns (Lovis 1977) and 13% of the taxa 
of Japanese pteridophytes with known mode of 
reproduction are reported to exhibit apogamous 
reproduction (Takamiya 1996). In apogamous 
species, all offspring from a parent are expected 
to be clonal with a small amount of genetic varia¬ 
tion. Interestingly, Darnaedi et ah (1990) report¬ 
ed that Dryopteris yakusilvicola Sa. Kurata, a 
triploid apogamous species of recent hybrid ori¬ 
gin, showed no allozyme variation within the 56 
individuals they examined. They considered D. 
yakusilvicola to have originated from a single hy¬ 
brid event between the tetraploid sexual species 
D. sparsa and the diploid sexual species D. sa- 
baei. 

This phenomenon, however, is not common. 
Despite the clonal nature of apogamous repro¬ 
duction, many apogamous species of ferns exhib¬ 


it extensive morphological and genetic variation 
and often form species complexes in which mor¬ 
phological species are difficult to distinguish be¬ 
cause of their continuous variation (Watano & 
Iwatsuki 1988, Suzuki & Iwatsuki 1990, Lin et ah 
1995, Grusz et ah 2009, Hori et ah 2014). This is 
because apogamous species of ferns can hybrid¬ 
ize with other sexual species to produce apoga- 
mously reproducible offspring (Lin et ah 1995, 
Ebihara et ah 2012, Lee & Park 2013, Hori et ah 
2014). 

Lee and Park (2013) recently attempted to as¬ 
sess reticulate evolution in subg. Erythrovariae 
sect. Variae sensu Fraser-Jenkins (1986; the 
Dryopteris varia complex) in Korea. They re¬ 
ported that D. sacrosancta (subg. Erythrovariae, 
sect. Variae) and D. chinensis (subg. Dryopteris, 
sect. Aemulae Fraser-Jenk.) shared the same chlo- 
roplast DNA (cpDNA) haplotype. Subsequently, 
Hori et ah (2014) reported triploid apogamous in¬ 
dividuals of the D. varia complex to have three 
PgiC sequences, which can be considered homeo- 
logs, on homeologous chromosomes, but not on 
homologous chromosomes. They reported that D. 
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chinensis and Dryopteris sacrosancta and D. ko- 
bayashii, which Serizawa (2009) separated from 
D. sacrosancta s.L, shared the same nuclear PgiC 
sequences. The data suggest that apogamous D. 
sacrosancta s.s. and D. kobayashii are hybrids 
between the D. varia complex and D. chinensis, 
which are distantly related based on the cpDNA 
tree (Ebihara 2011, Hori et al. 2014). 

Furthermore, Hori et al. (2014) reported that 
Dryopteris pacifica (subg. Erythrovariae, sect. 
Variae) and the D. erythrosora complex (subg. 
Erythrovariae, sect. Erythrovariae ) also shared 
the same PgiC sequences, suggesting that hybrid¬ 
ization with members of the D. erythrosora com¬ 
plex contributed to the diversification of the 
apogamous species of the D. varia complex. 
Thus, these results indicate that phylogenetically 
wider taxon sampling is necessary to understand 
reticulate evolution in apogamous ferns. 

Dryopteris hondoensis Koidz is a triploid 
apogamous species of sect. Erythrovariae, dis¬ 
tributed in China, Korea and Japan (Honshu, Shi¬ 
koku, Kyushu). It is distinguished from other 
members section Erythrovariae by its yellowish 
green, dull lamina and mostly plane scales on the 
rachis of the pinnae (Iwatsuki 1995). In this 
study, we determined D. hondoensis to be com¬ 
posed of nuclear DNA sequences from the D. 
erythrosora complex (D. erythrosora, D. caudip- 
inna and D. koidzumiana ) of sect. Erythrovariae 
and those from D. chinensis of sect. Aemulae. 
This is the first report of reticulate evolution be¬ 
tween subg. Erythrovariae sect. Erythrovariae 
and subg. Dryopteris sect. Aemulae. 

Materials and Methods 


Plant Materials 

Leaf samples were collected from six individ¬ 
uals of Dryopteris hondoensis and from possibly 
related sexual species in Dryopteris sections Ae¬ 
mulae, Erythrovariae and Variae, mainly from 
our field surveys in nine prefectures in Japan 
(Appendix). Their ploidy levels and reproductive 
modes are shown in Table 1. The reproductive 


mode of several individuals in the samples was 
estimated by counting the number of spores per 
sporangium (32, apogamous; 64, sexual) accord¬ 
ing to Manton (1950). All voucher specimens are 
deposited in MAK and/or TNS. 

Ploidy analyses 

Ploidy analyses followed the methods of Hori 
et al. (2014). To determine ploidy level, the DNA 
content (2C value) of each nucleus extracted from 
fresh pinnae was measured once per sample by 
flow cytometry using a Cyflow Ploidy Analyzer 
PA-II (Partec, Munster, Germany) and a Cystain 
UV Precise P kit (Partec). Approximately 100 
mm 2 of each pinna was torn into several pieces 
and finely chopped with a razor blade in 0.25 mL 
of nucleus extraction buffer from the kit. Then, 
0.8 mL of staining solution from the kit was add¬ 
ed to the chopped tissues. The crushed tissue and 
buffers were filtered through a 30-pm nylon mesh 
(Partec) before analysis. Approximately 25 mm 2 
of fresh leaf tissues of Nicotiana tabacum L. (2C 
value = 11.71 pg., Narayan 1987) was used as an 
internal standard. 

Molecular analyses 

For molecular analyses, a small amount of 
leaf tissue was dried in plastic bags, 20 cm * 10 
cm, using silica gel. Subsequently, total DNA was 
extracted from the dried leaves using cetyltri- 
methylammonium bromide solution according to 
the method of Doyle & Doyle (1987). 

The plastid gene rbcL was used in this study 
as the cpDNA marker. Polymerase chain reaction 
(PCR) amplification of rbcL was performed using 
the primers aF and cR of Hasebe et al. (1994). 
PCR entailed an initial denaturation step at 95°C 
for 10 min; followed by 35 cycles of denaturation, 
annealing, and elongation steps at 98°C for 10 s, 
50°C for 10 s, and 72°C for 7 s, respectively; and 
a final extension step at 72°C for 7 min. 

The nuclear gene PgiC, including exons 14- 
16 and introns 14-15, and GapCp (short) exons 
9-10 were chosen as the nuclear DNA (nrDNA) 
markers. The PgiC fragment was amplified using 
the primers 14F and 16R of Ishikawaeta/. (2002). 
The GapCp fragment was amplified using the 
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Table 1. The ploidy levels and reproductive modes of the materials species in this study. 


Species 

Ploidy 

Reproductive mode 

References 

Dryopteris hondoensis 

triploid 

apogamous 

Hirabayashi (1966, 1967), Mitui (1967, 1968) 

D. erythrosora 

triploid 

apogamous 

Hirabayashi (1966, 1967, 1974), Kurita (1961), Mitui (1968) 

D. caudipinna 

diploid 

sexual 

Hirabayashi (1970) 

D. koidzumiana 

diploid 

sexual 

Hirabayashi (1969), Mitui (1967, 1968) 

D. kinkiensis 

tetraploid 

sexual 

Hirabayashi (1967) 

D. chinensis 

diploid 

sexual 

Hori et al. (2015a) 


triploid 

apogamous 

Hirabayashi (1966, 1974), Kurita (1961), Mitui (1968) 


tetraploid 

apogamous 

Weng (1989), Nakato et al. (1995) 

D. protobissetiana 

diploid 

sexual 

Hori et al. (2015b) 

D. saxifraga 

diploid 

sexual 

Hirabayashi (1967), Lin et al. (1995), Mitui (1975) 


triploid 

apogamous 

Lee et al. (2006) 

D. sordidipes 

diploid 

sexual 

Hirabayashi (1969, 1974), Mitui (1967, 1968) 

D. varia 

diploid 

sexual 

Ebihara et al. (2014) 


diploid 

apogamous 

Hirabayashi (1966, 1967, 1974), Lin et al. (1995) 


triploid 

apogamous 

Hirabayashi (1970), Mitui (1966, 1968) 


tetraploid 

sexual 

Tsai & Shieh (1975, 1985) 

D. gymnophylla 

diploid 

sexual 

Kurita (1966), Hirabayashi (1969) 

D. sabaei 

diploid 

sexual 

Darnaedi & Iwatsuki (1987), Hirabayashi (1966, 1974), 
Kurita (1962), Mitui (1968) 


primers 132F (5'-GTGCTTCCGGAGTTA- 
AATGG-3') and 488R (5'-CAACATCATCTTC- 
GGTGTATCC-3'), which were newly developed 
in this study. PCR entailed an initial denaturation 
step at 95°C for 3 min; followed by 40 cycles of 
denaturation, annealing, and elongation at 98°C 
for 10 s, 58°C for 5 s, and 72°C for 5 s, respec¬ 
tively; and a final extension step at 72°C for 1 
min. PCR amplification was performed using 
PrimeSTAR Max DNA Polymerase (Takara, 
Kyoto, Japan) on a Model 9700 thermal cycler 
(Applied Biosystems, Foster City, CA, USA). 

The nucleotide sequences of rbcL were deter¬ 
mined by direct sequencing. For sequencing 
rbcL, aF, aR, cR (Hasebe et al. 1994), and D. pa- 
ci-bf (Hori et al. 2014) primers were used. 

PCR-single-strand conformation polymor¬ 
phism (SSCP) analysis was performed to exam¬ 
ine allelic variation at the nuclear marker level for 
each individual, following the method of Hori et 
al. (2014). Electrophoresis was performed using 
MDE gel solution (Lonza) with 2% glycerol at 
18°C for 16 h at 350 V for PgiC and 2% glycerol 
at 15°C for 9.5 h at 300 V for GapCp. 


To sequence the bands separated on the SSCP 
gels, pieces of the DNA bands in the dried gel 
were peeled off and incubated in 50 pL of TE 
buffer (10 mM Tris-HCl and 1 mM EDTA, pH 
8.0) at 37°C overnight. The supernatant solution 
was used as a template for further PCR amplifica¬ 
tion. Before sequencing, the PCR products were 
purified using lllustra ExoStar 1-Step (GE Health¬ 
care, WI, USA) and used as templates for direct 
sequencing. The reaction mixtures for sequenc¬ 
ing were prepared using a BigDye Terminator 
v.3.1 Cycle Sequencing Kit (Applied Biosystems). 
The reaction mixtures were analyzed using an 
AB1 3130 Genetic Analyzer (Applied Biosys¬ 
tems). 

All plant samples were classified based on 
their PCR-SSCP banding patterns. The genomic 
constitution of each band pattern was identified 
by determining the nucleotide sequences of each 
DNA band separated on the SSCP gel. 

Phylogenetic analyses 

For phylogenetic analyses, only one sequence 
representing each haplotype for cpDNA (rbcL) 
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and each allele for nrDNA ( GapCp , PgiC) was 
used in our datasets. The cpDNA and nrDNA se¬ 
quences were aligned using MUSCLE Citation 
(Edgar 2004) and analyzed separately by neigh¬ 
bor-joining (NJ), maximum parsimony (MP), or 
maximum likelihood (ML) analyses using 
MEGA version 6 (Tamura et al. 2013). The NJ 
tree was obtained using the p-distance method 
(Nei & Kumar 2000). The data are expressed as 
the number of base differences per site. All am¬ 
biguous positions were removed from each se¬ 
quence pair. The MP tree was obtained using the 
subtree-pruning-regrafting algorithm (Nei & Ku¬ 
mar 2000) at search level 1, in which the initial 
trees were obtained by the random addition of se¬ 
quences (10 replicates). In the ML analysis, the 
best-fitting model of sequence evolution for each 
DNA region was selected using MEGA version 6 
(Tamura et al. 2013). The GapCp tree was con¬ 
structed using the T92+G model, the PgiC tree 
was constructed using the HKY model, and the 
rbcL tree was constructed using the K2+G mod¬ 
el. The initial trees for the heuristic search were 
obtained automatically by applying neighbor-join 
and BioNJ algorithms to a matrix of pairwise dis¬ 
tances estimated using the maximum composite 
likelihood approach and then selecting the topol¬ 
ogy with superior log likelihood value. The indels 
were treated as missing characters in the MP and 
ML analyses. The bootstrap method with 1000 
replications was employed to estimate the confi¬ 
dence levels of monophyletic groups. Dryopteris 
sabaei (sect. Nephrocystis, Fraser-Jenk.), D. fra- 
grans (quoted from DDBJ) and two species of 
Polystichum (quoted from DDBJ) were used as 
outgroups. 

Morphological observation 

The morphology of Dryopteris hondoensis, 
D. chinensis, D. erythrosora, D. caudipinna and 
D. koidzumiana was compared, particularly fo¬ 
cusing on differences in the color and sheen of 
the lamina, shape of the lowest basiscopic pin¬ 
nule, and shape of the scales. 


Results 

Ploidy analyses 

The genome size of each species in this study 
was estimated as follows: Dryopteris koidzumi¬ 
ana (diploid sexual species; Hirabayashi 1969, 
Mitui 1967,1968), 17.14 pg (N= 1); D. caudipinna 
(diploid sexual species; Hirabayashi 1970), 18.04 
± 0.19 pg (N= 12); D. erythrosora (triploid apoga- 
mous species; Hirabayashi 1966, 1967, 1974, Ku- 
rita 1961, Mitui 1968), 23.48 ± 0.94 pg (N= 5); D. 
hondoensis (triploid apogamous species; Hira¬ 
bayashi 1966, 1967, Mitui 1967, 1968), 21.44 ± 
0.53 pg (N= 4); and D. kinkiensis (tetraploid sex¬ 
ual species; Hirabayashi 1967), 28.07 ± 0.92 pg 
(N = 3). The genome sizes of D. erythrosora and 
D. hondoensis (only the triploid apogamous cyto- 
type has been reported for this species in Japan, 
Table 1) were less than two-thirds of those of D. 
caudipinna and D. koidzumiana (diploid sexual 
species), but they were larger and smaller than 
those of D. kinkiensis (tetraploid sexual species). 
Those samples were therefore estimated to be 
triploid apogamous. 

Sequence variation of rbcL and molecular phylo¬ 
genetic trees based on variation 

In Dryopteris hondoensis and possibly related 
species, 12 types of rbcL sequences were found 
(Fig. 1, Appendix). The GenBank accession num¬ 
bers of the obtained sequences are shown in Ap¬ 
pendix. Among the 1205 sites, 105 (9%) were 
variable and 42 (3%) were parsimoniously infor¬ 
mative. NJ, MP, and ML analyses uncovered phy¬ 
logenetic trees with similar topology. The ML 
tree (highest log likelihood = -2451.5660) with 
bootstrap percentages (BPs) of NJ/MP/ML analy¬ 
ses are shown in Fig. 1. Dryopteris hondoensis, 
D. erythrosora, D. caudipinna, and D. koidzumi¬ 
ana displayed rbcL sequences belonging to the 
same clade. 

Sequence variation of PgiC and phylogenetic 
trees base on variation 

In each of the samples, several different se¬ 
quences of nuclear PgiC were detected by SSCP 
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1. ML tree (highest log likelihood = —2451.5660) based on sequence variation of plastid gene rbcL with BPs (>70) of NJ/ 
MP/ML analyses. 
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analyses. The number of different sequences was 
always the presumed ploidy level of the sample or 
less. Twenty seven different sequences of PgiC 
were identified in the samples of Dryopteris hon¬ 
doensis and in other species of Dryopteris ana¬ 
lyzed in the present study. The length of the se¬ 
quences varied from 607 to 661 bp. The data ma¬ 
trix for phylogenetic analyses included 672 char¬ 
acters after editing, of which 153 (23%) were 
variable and 68 (10%) were parsimoniously infor¬ 
mative. The ML tree (highest log likelihood = 
-2061.8095) with BPs of NJ/MP/ML analyses is 
shown in Fig. 2. Two monophyletic groups (A and 
B) were recognized, each of which was supported 
by a high BP. Group A contained the PgiC se¬ 
quences of D. erythrosora, D. caudipinna, and D. 
koidzumiana, whereas Group B contained se¬ 
quences of D. chinensis. Dryopteris hondoensis 
contained both group A and group B sequences. 
Some samples of D. hondoensis had two sequenc¬ 
es of group A and one sequence of group B; oth¬ 
ers had one sequence of group A and one se¬ 
quence of group B (Table 2). 

Sequence variation of GapCp and the phyloge¬ 
netic tree based on the variation 

In most of the samples, several different se¬ 
quences of nuclear GapCp were detected by 
SSCP analyses. In total, 21 different sequences 
were identified in the samples of Dryopteris hon¬ 
doensis and some Dryopteris species analyzed in 
the present study. The length of the sequences 
varied from 301 to 350 bp. The data matrix for 
phylogenetic analyses included 360 characters af¬ 
ter editing, of which 87 (24%) were variable and 
48 (13%) were parsimoniously informative. The 
ML tree (highest log likelihood = —1166.5740) 
with BPs of NJ/MP/ML analyses is shown in Fig. 
3. Four monophyletic groups (A, B, C, and D) 
were recognized, each of which was supported by 
a high BP. Group A contained the GapCp se¬ 
quences of D. chinensis, and Groups B-D con¬ 
tained those of D. erythrosora, D. caudipinna, 
and D. koidzumiana. Dryopteris hondoensis had 
three sequences, one each from groups A, B, and 
D (Table 2). D 2 and D 3 were distinguished by the 
positions of indels, but they are not distinguished 


on the molecular tree because indels were exclud¬ 
ed in the phylogenetic analyses. 

The estimated genotype of each sample is 
summarized in Appendix. It must be noted that 
our current method cannot distinguish differenc¬ 
es in gene dosage in polyploids; therefore, A^Q, 
A,C|, and A]C,C| were not distinguishable and 
were thus indicated as A^*. 

Morphological observation 

The morphology of Dryopteris hondoensis, 
D. chinensis, and species in the D. erythrosora 
complex (D. erythrosora, D. caudipinna and D. 
koidzumiana) is shown in Figs. 4-6 and Table 3. 
The adaxial lamina surface of D. hondoensis is 
yellowish green and dull, fresh green and dull in 
D. chinensis, and dark-green and shiny in the D. 
erythrosora complex (Fig. 4). The most basal ba- 
siscopic pinnule of the lowest pinnae of D. hon¬ 
doensis is slightly more shortened than the sec¬ 
ond-most basiscopic pinnule. Those of the D. 
erythrosora complex is clearly shortened than the 
second most basiscopic pinnule. Those of D. chi¬ 
nensis is more elongated than the second most ba¬ 
siscopic pinnule. (Fig. 5). Iwatsuki (1995) de¬ 
scribed the scales of D. hondoensis as being 
mostly plane, with bullate scales intermixed on 
the rachis of the pinna rachis. However, this is in¬ 
correct. These two types of scales are not inter¬ 
mixed on the pinna rachis. The scales of D. hon¬ 
doensis are bullate distally on the rachis of the 
pinnae (Fig. 6c), but plane on the basal portion of 
the pinnae (Fig. 6g). 

Discussion 

Dryopteris hondoensis has never been con¬ 
sidered to be closely related to D. chinensis. Fras- 
er-Jenkins (1986) included D. chinensis in 
Dryopteris subgen. Dryopteris sect. Aemulae, 
whereas he included D. caudipinna, D. erythro¬ 
sora, D. koidzumiana, and D. hondoensis in 
Dryopteris subgen. Erythrovariae sect. Erythr- 
ovariae. Thus, he (Fraser-Jenkins, 1986) did not 
discuss the similarity in morphological charac¬ 
teristics between these species. 

Our PgiC and GapCp analyses, however, sug- 
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Fig. 3. ML tree (highest log likelihood = —821.3017) based on sequence variation of nuclear gene GapCp with BPs (>70) of NJ/ 
MP/ML analyses is shown. 
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Table 2. Ploidy level, reproductive mode and constitution of nuclear maekers ( PgiC, GapCp). 


Species 

PgiC 

type 

GapCp 

type 

Ploidy level, reproductive 
mode 

number of 
samples 

locality 

D. caudipinna 1 

AgAg 

B,D, 

2x / sexual 

1 

Hyogo pref. 

D. caudipinna 2,3 

AfiAg 

b,d 2 

2x / sexual 

2 

Hyogo pref. 

D. caudipinna 4 

A4A4 

B,B, 

2x / sexual 

1 

Kanagawa pref. 

D. caudipinna 5 

A3A3 

B,B, 

2x / sexual 

3 

Kanagawa pref. 

D. caudipinna 6,7,12 

A 3 A<5 

B 2 B 2 

2x / sexual 

3 

Kanagawa pref. 

D. caudipinna 9 

AfiAg 

B,D, 

2x / sexual 

1 

Kanagawa pref. 

D. caudipinna 10 

AfiAg 

B,B, 

2x / sexual 

1 

Kanagawa pref. 

D. caudipinna 11 

A 3 A^ 

B,B, 

2x / sexual 

1 

Kanagawa pref. 

D. caudipinna 13 

A 3 A 8 

B,B, 

2x / sexual 

1 

Kanagawa pref. 

D. caudipinna 14,16 

A3A3 

BA 

2x / sexual 

2 

Kanagawa pref. 

D. koidzumiana 1 

A 8 A 8 

B,B, 

sexual 

1 

Kagoshima pref., Yakushima 

D. koidzumiana 2 

A7A9 

c,c, 

2x / sexual 

1 

Kagoshima pref., Yakushima 

D. koidzumiana 3 

AgAg 

B,C, 

sexual 

1 

Kagoshima pref., Yakushima 

D. koidzumiana 4 

AfiAg 

B,C, 

sexual 

1 

Kagoshima pref., Yakushima 

D. koidzumiana 5 

A9A9 

B,B, 

sexual 

1 

Kagoshima pref., Yakushima 

D. erythrosora 1 

As/Ag 

b,/d 2 

3x 

1 

Miyazaki pref. 

D. erythrosora 2 

AiA 5 A6 

b,d,d 3 

3x 

1 

Wakayama pref. 

D. erythrosora 3 

a,a s a 6 

B,/D, 

3x 

1 

Wakayama pref. 

D. erythrosora 4 

AjAA 

B 2 D]D 2 

3x 

1 

Wakayama pref. 

D. erythrosora 5 

A2A4A5 

b 2 d,d 2 

3x / apogamous 

1 

Tokyo-to. 

D. hondoensis 1,5,6 

a 2 a 5 b. 

a 2 b 2 d 2 

3x / apogamous 

3 

Tokyo-to., Wakayama pref. 

D. hondoensis 2,4 

A2A5B1 

A3B1D2 

3x 

2 

Tokyo-to. 

D. hondoensis 3 

A 2 /B, 

a 2 b 2 d 2 

3x 

1 

Miyazaki pref. 

D. chinensis 1 

B,B : 

A3A3 

2x / sexual (Hori et al. 2015) 

1 

Miyazaki pref. 

D. chinensis 4-6 

B 2 B 2 

A,A, 

2x / sexual (Hori et al. 2015) 

3 

Miyazaki pref. 

D. chinensis 8-10 

b 2 b 2 

A3A3 

2x / sexual (Hori et al. 2015) 

3 

Miyazaki pref. 

D. chinensis 1-3 

b,/b 2 

Aj 

apogamous 

2 

Tokyo-to., Wakayama pref. 


Table 3. Morphological differences between Dryopteris hondoensis and its related species. 



D. hondoensis 

D. chinensis 

D. caudipinna 

D. koidzumiana 

D. erythrosora 

Color of lamina 

yellowish 
green, dull 

fresh green, 
dull 

dark green, 
shiny 

dark green, 
shiny 

dark green, 
shiny 

The most basal basiscopic 
pinnule of the lowest 
pinna 

elongate a little 
than the 
second one 

elongate than 
the second one 

shorten than the 
second one 

shorten than the 
second one 

shorten than the 
second one 

Scales on distal positions 
of pinna rachis 

bullate 

flat 

bullate 

bullate 

bullate 

Scales on basal positions 
of pinna rachis 

flat 

flat 

bullate 

bullate 

bullate 
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Fig. 4. Color of pinnule of the second most basal pinna, (a: D. caudipinna, b: D. erythrosora, c: D. hondoensis , d: D. chinensis', 
scale = 1 cm) 



Fig. 5. Most basal basiscopic pinnule of lowest pinnae. (1: D. caudipinna , 2: D. erythrosora , 3: D. hondoensis , 4: D. chinensis ; 
scale = 1 cm). 
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Fig. 6. Scales on distal (a-d) and basal (e-h) positions of ra- 
chis of pinna, (a/e: D. caudipinna, b/f: D. erythrosora, 
c/g: D. hondoensis, d/h: D. chinensis ; scale = 1 mm) 


gested that Dryopteris hondoensis is of hybrid 
origin between Dryopteris. caudipinna, D. koid- 
zumiana, D. erythrosora and D. chinensis be¬ 
cause it contains DNA sequences matching the 
clades of those species in the molecular phyloge¬ 
netic trees (Figs. 2 & 3). Furthermore, D. hon¬ 
doensis possibly originated at least three times 
because D. hondoensis comprises three clones. 

The results of the maternally inherited rbcL 
suggests that the maternal parent of D. hondoen¬ 
sis was D. erythrosora (triploid apogamous), D. 
caudipinna and D. koidzumiana (diploid sexual 
species). Flori et al. (2014) suggested that diploid 
apogamous offspring originating from triploid 
apogamous individuals through unequal meiosis 
can be both maternal and paternal species in the 
reticulate evolution of the D. varia complex. 
Therefore, we have two hypotheses about the ori¬ 
gin of D. hondoensis. First, triploid apogamous 
D. chinensis (paternal) may result in unequal 
meiosis and then cross with diploid sexual D. 
caudipinna or D. koidzumiana (maternal). In this 
case, there may have been backcrossing to D. 
caudipinna, D. koidzumiana and D. erythrosora, 
because D. hondoensis shows two sequences be¬ 
longing to the D. erythrosora complex clades in 
the trees. Second, triploid apogamous D. erythro¬ 
sora (maternal) may have unequal meiosis and 
then cross with diploid sexual D. chinensis (pa¬ 
ternal). Additional nuclear markers are needed to 
determine the more plausible hypothesis. 


Dryopteris hondoensis and D. chinensis dis¬ 
play clear ecological and morphological differen¬ 
tiation; specifically, D. hondoensis is evergreen 
and has ovate to broadly ovate-subdeltoid lami¬ 
na, whereas D. chinensis is summer green and 
has pentagonal lamina. They have similar mor¬ 
phological characteristics, however, that coincide 
with the results from our DNA analyses that sug¬ 
gest a hybrid origin for D. hondoensis. Dryopter¬ 
is chinensis has fresh green lamina and lanceo¬ 
late scales on the stipes while D. hondoensis has 
yellowish green lamina and lanceolate scales on 
the stipes and basal portion of the rachis of the 
pinnae. Iwatsuki (1995) did not postulate a hybrid 
origin for D. hondoensis between the D. erythro¬ 
sora complex and D. chinensis, even though he 
listed several natural hybrids as being in Japan. 

In our study, we suggest the importance of 
phylogenetically wider sampling, including dis¬ 
tantly related species belonging to other sections, 
to clarify the origin of apogamous species of 
ferns and to resolve the reticulate relationships in 
Dryopteris. 

We are grateful to the following persons for their assis¬ 
tance in collecting plant materials: Dr. J. Yamashita 
(Okayama University); S. Kariyama (Kurashiki Museum 
of Natural History); and Mr. K. Akagi, Y. Inoue, T. Mina- 
mitani, K. Mizote, and K. Ohora (Nippon Fernist Club). 
This study was partly supported by a Grant-in-Aid for 
JSPS Fellows No. 26-1720 to K. H. and by a Grant-in-Aid 
for Scientific Research No. 25291089 to N. M. 
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Appendix. Voucher specimens examined in this study. Genotypes ( GapCp , PgiC ) identified by sequencing are in bold, other¬ 
wise, genotypes were deduced from comparison of band positions in SSCP gels. When ploidy is unknown, the genotype 
is in brackets. For samples with unknown genome dosage, unidentified genomes are in italic. *Quoted from Hori et al. 
(2014), **Quoted from Hori et al. (2015a), ***Quoted from Hori et al. (2015b). 


D. caudipinnaX 
D. caudipinnal 
D. caudipinndi 
D. caudipinnal 
D. caudipinnaS 
D. caudipinnab 
D. caudipinnal 
D. caudipinnaS 
D. caudipinnal 
D. caudipinnaXO 
D. caudipinnaX X 
D. caudipinnaX2 
D. caudipinna 13 
D. caudipinnalX 
D. caudipinnaX5 
D. caudipinna 16 
D. koidzumiana 1 
D. koidzumiana2 
D. koidzumiana3 
D. koidzumiana4 
D. koidzumianaS 
D. erythrosoraX 
D. erythrosora2 
D. erythrosora3 


Locality 

DNA 

amount 

Ploidy 1ml ^f uctive 

rbcL type 

PgiC type 

GapCp type 

Herbarium no. 

MAK TNS 

Japan, Hyogo pref., Himeji-shi 


sex 

AB908383* 

A 6 A„(LC093995) 

B|D 2 (LC094033,LC094035) 

411423 

Japan, Hyogo pref., Himeji-shi 


sex 

AB908383* 

A^Lawsoos) 

BjD 2 (LC094034,LC094035) 

411424 

Japan, Hyogo pref., Himeji-shi 


sex 

AB908382* 

A6A6(LC093995) 

B 2 D 2 (LC094034,LC094035) 

411425 

Japan, Kanagawa pref., Zushi-shi 


sex 

AB908384* 

A 4 A 4 (LC093994) 

B,B,(LC094033) 

1190554 

Japan, Kanagawa pref., Zushi-shi 

18.40 

2x sex 

AB908384* 

AjA,(LC093993) 

B,B 2 (LC094033,LC094034) 

432377 

Japan, Kanagawa pref., Zushi-shi 

17.95 

2x sex 

LC094043 

AjA^LCOOSOOSyLCOgsgOS) 

B 2 B 2 (LC094034) 

432378 

Japan, Kanagawa pref., Zushi-shi 

17.96 

2x sex 

AB908384* 

AjA^LCOOSOOSyLCOgsgOS) 

B 2 B 2 (LC094034) 

432379 

Japan, Kanagawa pref., Zushi-shi 

16.73 

2x sex 

LC094043 

A 3 Aj(LC093993) 

B 2 B 2 (LC094034) 

432382 

Japan, Kanagawa pref., Zushi-shi 

18.54 

2x sex 

AB908384* 

A 6 A s (LC093995,LC093996) 

B 2 B 2 (LC094034) 

432385 

Japan, Kanagawa pref., Zushi-shi 

18.54 

2x sex 

LC094043 

A 6 A 6 (LC093995) 

B|B,(LC094033) 

432386 

Japan, Kanagawa pref., Zushi-shi 

17.57 

2x sex 

AB908384* 

AjA 6 (LC093993,LC093995) 

B|B,(LC094033) 

432387 

Japan, Kanagawa pref., Zushi-shi 

17.57 

2x sex 

AB908384* 

AjA^LCOOSOOSyLCOgSOOS) 

B 2 B 2 (LC094034) 

432388 

Japan, Kanagawa pref., Zushi-shi 

17.84 

2x sex 

LC094043 

AjA g (LC03993,LC093996) 

B|B,(LC094033) 

432390 

Japan, Kanagawa pref., Zushi-shi 

18.10 

2x sex 

AB908384* 

A 3 Aj(LC093993) 

B|B,(LC094033) 

432392 

Japan, Kanagawa pref., Zushi-shi 

17.84 

2x sex 

AB908384* 

A 3 A 3 (LC093993) 

B 2 B 2 (LC094034) 

432397 

Japan, Kanagawa pref., Zushi-shi 

19.52 

2x sex 

AB908384* 

A 3 A 3 (LC093993) 

B,Bi(LC094033) 

432401 

Japan, Kagoshima pref., Yakushima Is. 


sex 

AB908381* 

A„A lt (LC093999) 

B|B,(LC094033) 

411078 

Japan, Kagoshima pref., Yakushima Is. 

17.14 

2x sex 

AB920329* 

A 7 A„(LC093998,LC094000) 

C,C,(LC094014) 

413260 

Japan, Kagoshima pref., Yakushima Is. 


sex 

AB908382* 

A S A,(LC093997,LC094000) 

B,C, (LC094014,LC094039) 

413262 

Japan, Kagoshima pref., Yakushima Is. 


sex 

AB920329* 

A 6 A„(LC093997) 

B,C, (LC094014,LC094040) 

413268 

Japan, Kagoshima pref., Yakushima Is. 


sex 

AB920329* 

A,A„(LC094000) 

B,B,(LC094039) 

413349 

Japan, Miyazaki pref., Nobeoka-shi 

23.42 

3x 

LC094041 

As/A 6 (LC094004,LC094005) 

B 2 /D 2 (LC094011 .LC094013) 

410287 

Japan, Wakayama pref., Shinguu-shi 

26.12 

3x 

LC094041 

A, A 5 A 6 (LC094001 ,LC094004,LC094005) 

B 2 D,D 3 (LC094009,LC094010,LC09401 

3) 

410977 

Japan, Wakayama pref., Shinguu-shi 

22.40 

3x 

LC094041 

A 1 A 5 A 6 (LC094001,LC094004,LC094005) 

B,/D 2 (LC094011 ,LC094912) 

410982 


D. erythrosoraX 
D. erythrosorab 
D. hondoensisX 


Japan, Wakayama pref., Nachikatsuura-cho 20.66 3x 

Japan, Tokyo-to, Hachioji-shi 24.80 3x 

Japan, Tokyo-to, Inagi-sbi 




LC094041 A, A 5 A 6 (LC094001 ,LC094004,LC094005) B 2 D,D 2 (LC094009,LC094010,LC09401 411063 

LC094041 A 2 A 4 A 6 (LC094002,LC094003,LC094005) B 2 D,D 2 <LC094009,LC094010,LC09401 ^ 

LC094042 A 2 A 5 B,(LC094006,LC094007,LC094008) A 2 B 2 D 2 (LC094018,LC094037,LC09403 406g33 


D. hondoensis2 


Japan, Tokyo-to, Inagi-sbi 


LC094042 A 2 A 5 B](LC094006,LC094007,LC094008) A 3 B,D 2 (LC094019,LC094036,LC09403 4102 4 g 


D. hondoensis3 
D. hondoensisX 
D. hoitdoensisS 
D. hondoensisb 
D. chinensisX 
D. chinensis2 
D. chinensisl) 
D. chinensisX 
D. chinensisS 
D. chinensisb 
D. chinensisl 
D. chinensisS 
D. chinensisl 
D. chinensisXO 
D. kinkiensis 
D. kinkiensis 
D. kinkiensis 
D. kinkiensis 
D. kinkiensis 


Japan, Miyazaki pref., Kitakata-machi 

22.91 

3x 


LC094042 

A 2 /B,(LC094006,LC094008) 

A,B,D 2 (LC094018,LC09403 7.LC094038 
) 

410346 

Japan, Tokyo-to, Inagi-shi 

20.49 

3x 


LC094042 

A 2 A 5 B,(LC094006,LC094007,LC094008) 

A 3 B, D,(LC094019,LC094036,LC094038 
) 

410957 

Japan, Tokyo-to, Inagi-shi 

21.57 

3x 

apo 

LC094042 

A 2 A 5 B,(LC094006,LC094007,LC094008) 

A 2 B 2 D 2 (LC094018,LC094037,LC09403 

8) 

410959 

Japan, Wakayama pref., Shinguu-shi 

20.82 

3x 


LC094042 

A 2 A 5 B,(LC094006,LC094007,LC094008) 

A 2 B 2 D 2 (LC094018,LC094037,LC09403 

8) 

410978 

Japan, Tokyo-to, Inagi-shi 




AB908373* 

B,/B 2 (AB908891*/AB908892*) 

yJ 2 (LC094017) 

405131 

Japan, Tokyo-to, Hachioji-shi 




AB908373* 

B,/B, (AB908891 */AB908892*) 

/1,(LC094017) 

405220 

Japan, Wakayama pref., Kozagawa-cho 




AB908373* 

B,/B, (AB908891*/AB908892*) 

yJ ; (LC094017) 

405262 

Japan, Miyazaki pref., Kitakata-machi 

13.01** 

2x** 

sex** 

AB908373 

B,B, (AB908892) 

A|A,(LC094017) 

410320 

Japan, Miyazaki pref., Hinokage-cho 

13.85** 

2x** 

sex** 

AB908373 

B,B, (AB908892) 

A|A,(LC094017) 

410321 

Japan, Miyazaki pref., Nobeoka-shi 

14.5** 

2x** 

sex** 

AB908373 

B,B, (AB908892) 

A|A,(LC094017) 

410342 

Japan, Miyazaki pref., Nobeoka-shi 

13.56** 

2x** 

sex** 

AB908373 

B1B2 (AB908891 /AB908892) 

A 3 A 3 (LC094020) 

410350 

Japan, Miyazaki pref., Nobeoka-shi 

13.22** 

2x** 

sex** 

AB908373 

B 2 B, (AB908892) 

A 3 A 3 (LC094020) 

410354 

Japan, Miyazaki pref., Nobeoka-shi 

13.22** 

2x** 

sex** 

AB908373* 

B 2 B,(AB908892*) 

A 3 A 3 (LC094020) 

412500 

Japan, Miyazaki pref., Hyuga-shi 

13.56** 

2x** 

sex** 

AB908373 

B 2 B 2 (AB908892) 

A 3 A 3 (LC094020) 

410356 

Japan, Shizuoka pref., Kawanehonn-cho 




LC094044 

LC093908,LC093910 

LC094015,LC094016 

406837 

Japan, Nara pref., Kashihara-shi 



sex 

LC094044 

LC093909,LC093910 

LC094015,LC094016 

409127 

Japan, Saga pref., Saga-shi 

29.86 

Xx 

sex 

LC094044 

LC093908,LC093910 

LC094015,LC094016 

410378 

Japan, Saga pref., Saga-shi 



sex 

LC094044 

LC093908,LC093910 

LC094015,LC094016 

410379 

Japan, Wakayama pref., Nachikatsuura-cho 

27.60 

Xx 

sex 

LC094044 

LC093908,LC093910 

LC094015,LC094016 

411046 


D. kinkiensis 
D. sordidipes 
D. sordidipes 
D. sordidipes 
D. sordidipes 
D. sordidipes 
D. sordidipes 
D. protobissetiana 
D. protobissetiana 
D. protobissetiana 
D. protobissetiana 
D. protobissetiana 
D. protobissetiana 
D. protobissetiana 
D. protobissetiana 
D. protobissetiana 
D. protobissetiana 
D. saxifraga 
D. saxifraga 
D. saxifraga 
D. saxifraga 
D. saxifraga 
D. saxifraga 
D. saxifraga 
D. saxifraga 
D. saxifraga 
D. saxifraga 
D. saxifraga 
D. saxifraga 
D. saxifraga 
D. saxifraga 
D. saxifraga 
D. saxifraga 
D. saxifraga 
D. saxifraga 
D. varia 
D. sabaei 
D. gymnophylla 


Japan, Wakayama pref., Nachikatsuura-cho 26.77 
Japan, Kagoshima pref., Yakushima Is. 

Japan, Kagoshima pref., Yakushima Is. 

Japan, Kagoshima pref., Yakushima Is. 

Japan, Kagoshima pref., Yakushima Is. 

Taiwan, Taipei City, Urai 
Taiwan, Taipei City, Urai 


Japan, Kagoshima pref., Yakushima Is. 14.63* * 

Japan, Kagoshima pref., Yakushima Is. 15.61*** 

Japan, Kagoshima pref., Yakushima Is. 16.58*** 

Japan, Kagoshima pref., Yakushima Is. 15.05*** 

Japan, Kagoshima pref., Yakushima Is. 

Japan, Kagoshima pref., Yakushima Is. 15.22* * 

Japan, Kagoshima pref., Yakushima Is. 16.10*** 


Japan, Kagoshima pref., Yakushima Is. 
Japan, Kagoshima pref., Yakushima Is. 
Japan, Kagoshima pref., Yakushima Is. 
Japan, Tokyo-to, Okutama-machi 
Japan, Tokyo-to, Okutama-machi 
Japan, Tokyo-to, Okutama-machi 
Japan, Tokyo-to, Okutama-machi 
Japan, Tokyo-to, Okutama-machi 
Japan, Tokyo-to, Okutama-machi 
Japan, Tokyo-to, Okutama-machi 
Japan, Tokyo-to, Okutama-machi 
Japan, Tokyo-to, Okutama-machi 
Japan, Tokyo-to, Okutama-machi 
Japan, Aichi pref., Shitara-cho 
Japan, Akita pref., Senboku-shi 
Japan, Akita pref., Senboku-shi 
Japan, Akita pref., Senboku-shi 
Japan, Akita pref., Senboku-shi 
Japan, Akita pref., Senboku-shi 
Japan, Akita pref., Senboku-shi 
Japan, Akita pref., Senboku-shi 
Taiwan, New Taipei City, Pingling 
Japan, Tokyo-to, Inagi-shi 
Japan, Tokyo-to, Hachioji-shi 


Xx 


2x*** 

2x*** 

2x*** 

2x*** 

2x*** 
2x *** 


LC094044 LC093908,LC093910 

AB908372* AB908894* 

AB908372 LC093900,LC093901 

AB908372 LC093904,LC093905 

AB908372 AB908894 

AB908372 AB908894 

AB908372 LC094023,LC094024 

AB9I4474* AB914476* 

AB9I4474* AB914475* 

AB914474* AB914475* 

AB914474* AB914475* 

AB914474* AB914475* 

AB914474* AB914475* 

AB914474* AB914476* 

AB914474 AB914475 
AB914474 LC093906 

AB914474* AB914475* 

AB908375* AB908878* 

AB908375* AB908878* 

AB908375 AB908878 

AB908375 AB908878 

AB908375 AB908878 

AB908375 AB908878 

AB908375 AB908878 

AB908375 AB908878 

AB908375 AB908878 

AB908375 AB908878 

AB908375 AB908878 

AB908375 AB908878 

AB908375 AB908878 

AB908375 AB908878 

AB908375 AB908878 

AB908375 AB908878 

AB908375 AB908878 

AB908375 AB908878 

AB908376* AB908876* 

LC095874 LC093907 

AB908370* AB908893* 


LC094015,LC094016 411052 


LC094025 

LC094026 

LC094023.LC094026 

LC094023 

LC094025 

LC094024,LC094026 

LC094028,LC094029 

LC094028 
LC094028 
LC094028 
LC094028 
LC094028 
LC094028 
LC094028 
LC094028 
LC094028 
LC094030 
LC094030 
LC094031 
LC094030 
LC094030 
LC094031 
LC094031 
LC094030 
LC094030 
LC094030 
LC094030 
LC094030 
LC094030 
LC094031 
LC094030 
LC094030 
LC094030 
LC094030 
LC094027 

LC094021 ,LC094022 
LC094032 


407057 

410918 

410919 

410920 

423478 

423482 

410906 

410907 

410913 

410914 

410915 

410916 

410917 

425899 

417179 

405345 

405346 

405347 

405348 

405349 

411382 

411383 

411385 

411386 

411387 

411406 

423578 

423579 

423584 

423586 

423587 

423588 

423590 

406981 

405219 











